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Abstract: We outline a new formalism for dipole-type parton showers which main- 
tain exact energy-momentum conservation at each step of the evolution. Particular 
emphasis is put on the coherence properties, the level at which recoil effects do enter 
and the role of transverse momentum generation from initial state radiation. The 
formulated algorithm is shown to correctly incorporate coherence for soft gluon ra- 
diation. Furthermore, it is well suited for easing matching to next-to-leading order 
calculations. 
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1. Introduction 

Parton shower simulation programs like Pythia and HERWIG [1,2] have been the 
workhorses for high energy physics experiments for a long time. With the advent of 
the Large Hadron Collider (LHC) at CERN many new developments have been made 
in order to refine the existing programs and to extend their applicability. Currently, 
there are a few major multi purpose event generators in use. The FORTRAN gener- 
ators have been completely rewritten and extended as PythiaS and Herwig++ [3-5] 
and the program Sherpa has been established [6,7]. All these new generators come 
with a great effort of new developments. 

The classical parton showers with 1 — 2 like branchings themselves have been 
rewritten and reformulated in order to take into account mass effects and improve 
the phase space coverage [8] or to incorporate intertwined interactions with partons 
from additional hard interactions in the underlying event [9]. Other developments 
try to incorporate quantum interference effects or subleading terms in the expansion 
in the number of colours [10-13]. Besides many other efforts to improve the event 



simulation, e.g. in the area of the underlying event [14-19], one of the most important 
questions that has been addressed was for the matching with higher order matrix 
elements or the merging of parton showers with multiple hard emissions at the tree 
level. 

The latter was formulated conceptually by Catani, Krauss, Kuhn and Webber 
(CKKW) [20,21] and, in a variant, by Lonnblad [22] and implemented in various 
parton shower programs [23-27]. An alternative approach to matching has been 
advocated by Mangano [28, 29] and is known as MLM matching. [29] gives a com- 
prehensive overview of current implementations and shows detailed comparisons. 
Recent improvements address previous deficiencies of the CKKW algorithm in filling 
the phase space [30,31] that had been discussed extensively in [32]. Some problems 
have been overcome with the help of so-called truncated showers that were intro- 
duced earlier in the context of matching with next-to-leading order (NLO) matrix 
elements [33] . Merging with NLO matrix elements has been recently studied in [34] . 

First attempts to match parton showers and matrix elements at NLO have been 
made in a phase space slicing approach [35-37] that suffered from some (numerically 
small) systematic inconsistencies. The more general and systematic approach has 
been MCCilNLO [38,39] that has since been extended to include many processes and 
all possible colour structures [40-42]. As the MC@NLO is closely tied to a partic- 
ular parton shower algorithm, the subtraction terms that guarantee the consistent 
matching with the NLO matrix element contributions have to be calculated once for 
a specific parton shower program. After all initial efforts have been made to match 
with HERWIG, later also several processes have been matched to Herwig-I— |- [43-45]. 

While the MC@NLO approach is very successful and was developed for many 
processes it may be considered to be tied too strongly to the underlying parton 
shower and to suffer from negative weighted events, which, in practice, never pose 
a real problem. An alternative approach, now known as POWHEG [33,46-48], has 
been formulated by Nason [33]. Here, the matching formalism is based on a mod- 
ified Sudakov form factor that contains the real emission matrix element. This 
guarantees that the first emission of the parton shower is the hardest one as well 
and can, therefore, be described by the full matrix element. In principle, this ap- 
proach is closely related to the so called hard matrix element corrections [49,50]. 
This approach has been further developed into a systematic matching scheme and 
apphed to several processes [51,52] and was also widely used by other groups, e.g. 
the Herwig++ collaboration [45,53-56]. 

The subtraction based original MC@NLO approach is following the requirement 
that the underlying parton shower algorithm must not be modified. Already shortly 
after the publication of the MC@NLO approach it was, however, noted that the 
intrinsic subtraction in this scheme could be simplified tremendously if the parton 
shower would follow closely the subtraction terms that are used to regularize the 
soft and coUinear divergences in the NLO calculation to be matched. This was 
highlighted in [10], where also more formal developments towards new partons were 



carried out. Therefore, some groups started to write new parton showers based on 
subtraction terms. The groups [57, 58] used the Catani-Seymour (CS) [59] terms, 
while [60] based a shower on Lund dipoles and [61] on antenna subtraction terms. 

In this paper we take the same viewpoint and present new theoretical devel- 
opment towards a parton shower, based on CS subtraction kernels as well. After 
the question for soft coherence effects has been raised by [62] in the context of a 
simplified toy model, we would like to address this question for a full CS like par- 
ton shower. The question of collinear radiation and the implication of the DGLAP 
evolution was discussed in [63,64], not the coherence properties, though. In this 
paper, we formulate a parton shower, based on CS subtraction kernels, in a detail 
that specifies a full implementation. We show that, with the right choice of evolution 
variable and initial conditions it is indeed possible to find the correct soft anomalous 
dimensions with a parton shower, based on CS dipoles. Therefore, it is possible to 
implement such a shower while incorporating soft colour coherence effect in a way 
that has always been a vital ingredient of the HERWIG and Herwig-I— |- programs. 



2. Local Recoils, Form Factors and Coherence 

We consider a single parton emission off a pair of partons with momenta Pij and p^. 
The probability for this emission is taken to be the sum of two splitting functions, 
each associated with one leg. Using DGLAP splitting kernels and the Sudakov 
decomposition for the splitting pij — > qi, qj, 

Qi = zpij + n + k± , (2.1) 

2pij ■ n z 

■^. = (i-^)p. + ap,.f(i-.) "-*^- 

where — —p\ and ■ Pij — k± ■ n — constitutes the usual collinear approxi- 
mation, which may be extended to the quasi-coUinear approximation for emissions 
off massive partons, [8]. The light-like vector n defines the collinear direction, and 
therefore is used as the gauge vector in a light-cone gauge when deriving the coUinear- 
singular behaviour of QCD matrix elements, n needs to be chosen along the colour 
connected partner p^, the so-called physical gauge, in which interference diagrams 
are coUinearly subleading such that the unregularized splitting kernels are given by 
cut self-energy diagrams only.^ 

Note that, within this parametrization, the DGLAP splitting kernels are func- 
tions of 

n ■ qi 



n ■ Pij 



(2.3) 



^We note that this is a gauge choice for each singular Hmit of interest. The definition of 'colour- 
connected' here applies in the large-A^c limit but may be generalized by including the full colour 
correlations present at finite Nc- 



Indeed, there is not a single choice of hght-cone gauge, but rather a class of gauge 
choices which are connected by rescaling the gauge vector n {i.e. longitudinal boosts 
along the coUinear direction), for which the splitting kernels are left invariant. 

We are interested in extending this picture such as not to perform an approx- 
imation in the choice of kinematics, thereby introducing exact energy-momentum 
conservation within the splitting Pij,Pk — > Qi, qj, Qk- The choice of the recoil strategy 
is not unique. However, choosing a spectator to absorb the longitudinal recoil of the 
splitting, 

n^Pk = ( 1 - ^^T- r ) Pk (2.4) 

V 2pij ■ Pk z{l - z) J 

is the only choice compatible with the remaining gauge degrees of freedom in the 
functional form of the sphtting kernels. As we shall also see, this is the only choice 
which guarantees that the splitting functions in a physical gauge do reproduce the 
correct soft behaviour. 

2.1 DGLAP Kernels, 'Soft Correctness' and Angular Ordering 

As we are primarily interested in soft gluon radiation, we neglect gluon splittings 
into quark-antiquark pairs in this section. 

For final state radiation the spin-averaged DGLAP kernels are given by 

P,A^) = + (1 - ^)) , P,,{^) = 2C^ + + ^(1 _ , 

(2.5) 

such that matrix elements squared, summed over all coUinear configurations factorize 
as 

n+l {M{qi, qn+i)\M{qi, ?n+l))n+l 
™ 4 

,Pij,...,qn+i)\M{qi,...,Pij,...,qn+i))n ■ (2.6) 

i=l j^i ' 

Note that in writing this expression, we do need to include a symmetry factor of 1/2 
along with the gluon splitting function. 

As each amplitude \M) is a colour singlet, i.e. 

n n 

Et' + EEt^-T,- = (2.7) 

i=l 1=1 j^i 

we may rewrite coUinear factorization within the choice of the physical gauge for 
single coUinear configurations as 

{M{qi, qn+i)\M{qi, gn+l))n+l 

n 4,-jioi 

-Pij{z)\n=pk n{M{qi, ...,Pij, ...,qn)\C^ij)k\M(qi, ...,Pij, ...,qn))n , (2-8) 

i=l j,k^i ■ 



where 



T, ■ T, 



(2.9) 



is the colour correlation operator as introduced in [59] . 
Within this framework, we have that 

qi-Pk 11- 



Qi ■ Qj 1 



n=Pk 



Qi ■ qj Qj ■ Pk 



Qi • Qj 



Qj • Pk 



n=Pk 



Qj -QiQi- Pk 



(2.10) 



such that the single splitting function Pij{z)\n=pi_ constitutes the complete, correct 
soft behaviour for the dipole i,k. Note that the eikonal parts - as well as any 
other part of a splitting function - is invariant under rescaling of the spectator 
momentum pk, which is an even stronger motivation to use the longitudinal recoil 
strategy defined above. 

This will also be a necessary requirement when trying to remove what we call 
'soft double counting'. As we will show now, this is closely related to the coherence 
properties and logarithmic accuracy of a particular shower setup. To be precise, 
we consider the form factor Ajfe((5^,/U^) associated to a final-final dipole i.k when 
evolving from a hard scale to a soft scale /j?. Regarding the leading- (dou- 
ble) and next-to-leading (single) logarithmic contributions, a^L^" and q;"L^"~^ with 
L — ln((5^//i^) the correct behaviour can be obtained from the coherent branching 
formalism [65], reproducing the results of soft gluon resummation, [66], by consid- 
ering the leading behaviour of the 2;-integrated splitting kernel for /i^ <^ p\ <^ Q^. 
The resulting form factor reads 



-lnA,fc(g2,/.2)= / 

J u 



dpi asip] 



Pi 



2tt 



{T,{pI,Q')+T,{pI,Q')) , (2.11) 



where the Sudakov anomalous dimensions T^^pl^Q^) are given by 



^,{pI,Q') 
r.(pi,Q') 



pi 

pj_ 6 



(2.12) 
(2.13) 



receiving contributions both at the LL level from soft collinear, at the NLL level 
from hard collinear radiation. Note that the latter, i.e. the non-logarithmic terms 
in r are determined by the average of the soft-suppressed, 2;-regular terms of the 
splitting functions. 



2.2 Recoils and Soft Coherence 

We now want to include the effects of a finite recoil. Within the minimal recoil 
strategy outlined above this only affects the phase space measure. 



dpi 

pi 



■dz 



dp]_ 



d^ 1 - A 



2 

PI 



z{l - z)Sik 



Sik = • Pk 



(2.14) 



where we introduced A ^ 1 to explicitly keep track of these effects. Choosing a 
phase space region related to an ordering in virtuality or transverse momentum, 



V < < Q' , (2.15) 

z[l - z) 



we find 



Q' 3 ^ 
PI 2 



T:{pl,Q') = C^{2\n^-^--2X^] , (2.16) 



r^{pl,Q') = CA(2ln^-^-^-2X^) . (2.17) 

\ P± C) Sik J 

Note that here, the recoil effects enter at the level of next-to-leading logarithms and 
the coefficient of the leading logarithms turns out to be twice the correct result. 
The latter observation has been noted since long [65]. From this example it is very 
clear that the simple fact that the DGLAP splitting functions reproduce the correct 
soft behaviour is not enough for the correct soft anomalous dimension. The wrong 
coefficient of the leading logarithmic contributions may be attributed to a double 
counting of soft emissions, originating from the fact that the above chosen phase 
space region does introduce an overlap of the phase space available for emissions off 
either parton of the dipole. 

Choosing angular ordering in the variable q by disentangling soft and collinear 
limits^, and imposing phase space constraints through a cutoff on the transverse 
momentum in the soft limit (s), 

e = /^^ < , (1 - z)Y e < Q' (2.18) 

z-^il — zV 



we recover the correct anomalous dimensions (|2.12| , p.l3|) with recoil effects entering 
beyond NLL, 

rr (Pi. Q')-C. (m I - I) . (l - 2Ag) + O (I) . (2.9) 

rf (Pi. Q^) ^ (.„ fyi)^ 2C.| (l - Ag) + O (I) . (2.20) 

the subleading terms giving rise to power corrections in the form factor exponent. 

Apart from the recoil effects, this result has a straightforward explanation: The 
phase space region chosen for the angular ordered evolution provides approximately 
disjoint regions (exactly disjoint in the case of Herwig++) for emissions off either 
leg of the dipole, thereby removing the soft double counting observed earlier. Note 
that this observation would then in principle allow to include local recoils within the 
angular ordered DGLAP evolution. 



^Note that in the soft hmit(s), z ^ t , 1 — p\ scales as e^. 



2.3 Catani-Seymour Kernels and New Formalism 

As outlined in the previous sections, taking a minimal choice to treat recoils yields 
a dipole-type picture. Within such a cascade it is however difficult to maintain 
the strong angular ordering, which is tied to the 1 — > 2 nature of independent jet 
evolution. 

Choosing the phase space to be restricted by a cutoff on the transverse mo- 
mentum, thereby assuming an ordering in (or virtuality) is a much more natural 
picture to consider for a dipole-type evolution. In addition, this also removes com- 
plications when implementing matrix element corrections, either stand alone or for 
the purpose of POWHEG-type NLO matching as the first emission off a dipole then 
is indeed the hardest emission. 

To cure the problem of soft double counting generated by this evolution, one 
may modify the splitting functions and 'continue' them over the whole available 
phase space in such a way, that the soft-singular pieces reproduce the correct soft 
behaviour when adding both modified splitting functions. 

More precisely, for each leg i we replace the eikonal part by the radiation pattern 
associated with coUinear emissions of pi 

Pi ■ Pj _^ Pi • Pj /2 21) 

Pi-qpj-q Pi- q (Pi+Pj) ■ q 

while keeping the collincar parts exactly. Note that this minimal construction, which 
does not modify the singular properties following from QCD, is nothing but the 
construction prescription for the subtraction kernels introduced in [59] . This picture 
of local recoils using a single spectator parton is ideally supplemented with exact 
factorization of the phase space considering no kinematic approximation. One choice, 
which so far has been implemented [57, 58] is to invert the kinematic mappings as 
derived in [59]. 

For initial state radiation, taking the Catani-Seymour factorization literally does 
have shortcomings. Most prominently, the choice of keeping the initial state emitter's 
momentum collincar to the one before emission leads for example to the fact that 
a final state singlet as in Drell-Yan Icpton pair production, does receive a non- 
vanishing transverse momentum from the very first shower emission only. Further, 
an initial-initial system emitting a parton left the spectator parton unchanged, which 
might not be sufficient for the description of the transverse momentum spectrum of 
the whole final state. The aim of this work is to provide a formalism, which does 
overcome these problems. Further, we are interested in the logarithmic accuracy 
and ordering of soft gluon radiation in our setup reflecting coherence properties. 

Starting from the final-state parametrization given above, the outline of our 
formalism is as follows: We obtain a parametrization of the kinematics for initial 
state emitters and/or spectators by considering the physical splitting processes while 
maintaining exact energy-momentum conservation locally to each branching, i.e. 
involving the emitter-emission system and a single spectator only. The spectator 



is restricted to take the longitudinal recoil of the splitting only. For initial state 
radiation we do allow each initial state emission to generate transverse momentum of 
the emitting incoming parton in a backward evolution. This transverse momentum is 
then migrated to the complete final state system by realigning the incoming partons 
to the beam axes at the end of the evolution. 

For final-final dipoles, we find that the anomalous dimensions take the correct 
form apart from the fact, that the dependence on the arbitrary hard scale is 
being replaced by the dipole's invariant mass Sjfe, 

rf(pi.-)^C.(ln|-2)-C..Aj= + o(|) . (2.22) 



rf (pi, .) = C J In 5 - ii - C^ttA^ + O Mi ) , (2.23) 



^0 



itol^l . (2-24) 



with recoil effects entering beyond NLL. 

We note that, in case of DGLAP kernels, the correct coefficient of the leading 
logarithmic contributions to the anomalous dimension is governed by the choice of 
boundaries on the momentum fraction for a given (but arbitrary) hard scale Q^, 

"^-^ d^ _ 1 /g2~ 
with K = pi/Q^. 

The above findings for the anomalous dimension can essentially be traced back 
to the fact that the transition from a DGLAP kernel possessing a soft singularity 
~ 1/(1 — 2;) to the appropriate Catani-Seymour kernel (while keeping track of all 
recoil effects, i.e. considering the soft limit at fixed p^) is the simple replacement 

where Kjfc = p\/ Sik- Here, the first factor is the effect of the finite recoil stemming 
from the exact factorization of the phase space measure. 

Within the variables to be outlined in detail in the next section, we find that 
this pattern generalizes to the cases of initial state emitter or spectator partons, up 
to a sign on the recoil term owing to timelike or spacelike virtualities of the emitter 
or whether the relevant dipole scale is a spacelike momentum transfer or invariant 
mass. 

Choosing the z boundaries (in the approximation considered above) to be given 

by 

2 

z<l-^^l-K (2.26) 

it is evident that the recoil contribution only gives rise to power corrections, while 
the logarithmic contribution is given by 



2/,« 



^ In ( ^ ) + power corrections , (2.27) 



2 J«2 + 2 \p]_ 



thereby reproducing the correct coefficient up to the disappearance of the arbitrary 
hard scaie, an immediate consequence of the screening of the soft singularity at fixed 
transverse momentum. 



2.4 Structure of the Evolution 

For final state radiation with final state spectator, our findings of the previous section 
immediately signal a choice of the hard scale for a single dipole originating from a 
hard process. Choosing an arbitrary hard scale ^ Sik will immediately result 
in the appearance of spurious logarithmic contributions when performing the p\ 
integration. 

For example, at fixed the leading logarithmic contributions for a dipole i, k, 
with Casimir operators Ci^k associated to the partons, take the form 

- In ^ £(C.' + CI) m g) In (-iL) + NLL (2.28, 



instead of the expected result 



■ In A,, = ^{Cf + CD In^ ) + NLL , (2.29) 



the mismatch being manifest as an ambiguity at the level of next-to-leading loga- 
rithms. We are therefore lead to the choice — Sik, i.e. the hard scale associated 
to a dipole is the respective invariant mass. 

For initial state emitter or spectator partons, we assume that this generalizes to 
choosing the hard scale in such a way as to fill the complete phase space, modulo 
the infrared cutoff. 

Note that this choice does not determine the ordering per se, but only the choice 
of hard scale and the shape of the phase space restriction when evolving between 
two scales. The ordering is to be chosen in such a way, that the leading effects of 
multiple emissions off each leg of the dipole do exponentiate. Due to the structure 
of the splitting kernels given above and the additional complications from all finite 
recoil effects the explicit exponentiation is beyond the scope of this paper. 

Having however observed that we can reproduce the correct Sudakov anomalous 
dimension, while avoiding soft double counting we additionally note that within the 
variables chosen 

pI = 2^'-^^-^^ (2.30) 

Pi -Pk 

for emission of a gluon of momentum g off a dipole {i,k). Ordering emissions in this 
variable therefore corresponds to an ordering reproducing the most probable history 
of multiple gluon emission according to the eikonal approximation in the limit of 
soft gluons strongly ordered in energy. 

We therefore conclude that branchings within the physical kinematics outlined 
above and based on the corresponding CS dipole splitting functions allow us to 



construct a parton shower that has the right coherence properties. The final state 
emissions should in this case be taken as outlined above, i.e. with the hard scale 
of a single cascade chosen to be the dipole invariant mass and the evolution should 
be strictly ordered in transverse momentum. However, in the naive adoption of the 
CS picture to a parton shower not every initial state emission would contribute to 
the final state transverse momentum. We will formulate a more suitable approach 
below. 



3. Kinematics, Phase space and Splitting Probabilities 

The purpose of this section is to provide our new formalism in all details, particularly 
the kinematic parametrization, phase space convolution properties and boundaries, 
and the splitting probabilities entering the evolution. The explicit expressions for 
the splitting kernels can be inferred from [59], as we express the variables chosen 
there in terms of the physical variables chosen for the shower evolution. 

3.1 Final State Radiation 
3.1.1 Final State Spectator 

Final state radiation with a final state spectator docs represent the generic version 
of the splitting kinematics chosen here. For a splitting {pi,Pj) — > {qi, q, qj) we choose 
the standard Sudakov decomposition 

Qi = zpi H —pj + k± (3.1) 

ZSjj 

q^{l- z)pi + Pj - (3.2) 

= " z{l -".). J ' ^'-'^ 

where — — 2pj • pj = —Sij and k± ■ Pij — 0. The transverse momentum is defined 
in the dipole's rest frame to be purely spacelike, 

PiJ = (^-^'±P^ ' k± = {0,p±), P-P± = 0. (3.4) 

Note that this does preserve the momentum of the emitting system, qi + q + qj — 
Pi +Pj- The parametrization gives rise to the phase space factorization [59] 

1 d0 dz f P^ \ 

d03(g.g,g,|g) = d03(p.P.|g)^2^dpi^^ [l- ^^^_\^^J (3-5) 

Note that, in the coUinear limit, this is the massless version of the kinematics as 
chosen in [8]. It further constitutes the inversion of the 'tilde'-mapping, where the 
variables y and z chosen in [59] are given by 
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Figure 1: Allowed phase space regions for emissions from a final- final dipole expressed in 
the Dalitz variables Xk = 2Q-pk/Q'^ for a dipole of mass Sij = 100 GeV and infrared cutoff 
^ = f> GeV. The shaded region is accessible for emissions off the parton i, whereas the 
area enclosed by the solid line is accessible for emissions off parton j. The area enclosed 
by the dotted line is an example of the phase space excluded when starting at a scale lower 
than Sij. Note that the infrared cutoff is exaggerated for illustrative purposes only. In 
practice, almost the whole physical phase space will be available. 

The allowed phase space region is obtained by considering the limits on the 
emitter's virtuality, 

.2 ^ Pl ^ n,2 



such that 



4/^' < ;^fr7) < ^-x = (3.7) 




Averaging over azimuth, the final-final splitting kernels take the form 

^snp\.^y) (3.9) 

such that the splitting probability is 

dP^|Wpi.--)>(l-j^)fd.. (3.10) 

Note that, comparing to the coUinear limit, the effect of finite recoils is to act as 
a damping factor for large-angle hard emissions, provided that y < \ which is a 
consequence of the phase space boundary. 



3.1.2 Initial State Spectator 

For an initial state spectator we consider the crossing qj — > —qa, Pj — —Pa, such 
that 

Qi = zpi + ^Pa + k± (3.11) 

2 

g = (1 - z)pi + Pa - (3.12) 

g„= fl + )pa. (3.13) 

V z{l-z)SiaJ 

Note that exact momentum conservation is trivially implemented by just the fact 
that the parametrization for a final state spectator does respect this constraint. The 
transverse momentum is defined be purely spacelike in a frame where, 



Pi,a = (^^, ±PJ , = (0, p±) , p • p± = . (3.14) 

The phase space measure then obeys the convolution 

d(p2{(li:(l\Q;Pa:qa:Xa) = d(p( {pi\Q; Pa, Pa, XXa)^:^^^^^^—^dp'i , (3.15) 

where 

J- Pa ' Qi /o 1 /-•\ 

X — ^2 ) z — (3.lDj 

and it is straightforward to verify that this indeed gives rise to the phase space 
convolution as given in [59]. Including the parton distributions and the kinematic 
factor of the partonic flux, the relevant measure is 



Ma ■ n 



d02 {Qu Q\Q; Pa, Qa, Xa)dXa = 



0{X - Xa) — —-^-7Z rdp_L d0i {pi\Q;Pa,Pa,Xa)dXa . (3.17) 



fa{xa) 27r 167r2 z{l -z) J ^Pa-n 

The phase space hmits can be obtained as for the final state case, 

2 

V < < Ql^ , (3.18) 

where, owing to x > Xa, the hard scale of a dipole is now given by 

Vmax — ^ia " • (3.19) 

Xa 



Averaging over azimuth, the final-initial splitting kernels take the form 

^^"•^^npi,-)) (3.20) 



2qi-qx 



such that the sphtting probabihty is 

= SWx, z))lj^9{x - x^)^dz . (3.21) 

Note that the finite recoil enters only in the PDF ratio, reproducing the correct 
coUinear limit when x ^ 1. Once again, the efi^ect of the finite recoil is a damping 
of hard emissions for x Xa- 

3.2 Initial State Radiation 

A construction of initial state radiation by just crossing prescriptions is not obvious 
owing to the fact that the shower evolution is formulated as a backward evolution. 

The physical variables thus need to be defined from the physical forward kine- 
matics. For the physical emission process Qa — > Pa, Q the relevant Sudakov decom- 
position for the emission momentum q is 

^forward = (1 " ^)9a + Tj 7Z - k± , (3.22) 

Zn ■ qa[l - z) 

where n is the backward lightcone direction defining the coUinear direction, i. e. the 
final or initial state spectator's momentum. 

The parametrization above is most conveniently inverted to backward evolution 
Pa qa, q by considering the process in a frame where qa = Pa/x, giving rise to 

(1 — z) 1 

^backward = Pa + 7) _ 1=^^ ■ (3.23) 

X Zn • Pa\L Z) yX 



We therefore define the Lorentz invariant physical variables to be given by 

P^ 

X qa-q = X n ■ q = {1 - z)n ■ Pa . (3.24) 

1 — z 

The parametrization keeping the emitter aligned with the beam axis can then be 
related to a parametrization where the initial state parton after (backward evolution) 
emission does acquire a finite transverse momentum while keeping the spectator after 
emission aligned with the one before emission. 

It is this type of splitting kinematics which allows any emission off an initial 
state parton to contribute transverse momentum to a final state system after having 
applied a proper realignment boost once the parton shower evolution has terminated. 
Ideally, this final boost should not be related to the parametrization chosen but being 
determined in a process dependent way such as to leave the interesting kinematic 
quantities of the hard process invariant. 

3.2.1 Final State Spectator 

For initial state emissions with final state spectator, Pa,Pj — > qa,Q,Qj, using the 
variables introduced in [59], 

x= , ""-"^ «= , , (3.25) 

[Pa - Pj) ■ qa [Pa " Pj) ' ?a 



we use the parametrization 

1 —U U 1 — X 1 , /o r^^X 

Qa = Pa H Pj H k± (3.26) 

X ~ U X X — u u — X 

1 — X u 1 — u 1, ,_ 
q = Pa H Pj H k± (3.27) 

X — U X X — u u — X 

<l3= {^-1)P^^ (3-28) 

which does preserve the momentum transfer, q + qj — Qa — Pj — Pa- The transverse 
momentum obeys 

1 — X 

k\ = -u{l - u) Saj Saj = Ipa ' Pj ■ (3.29) 

Considering the colhnear hmit -u 0, it is evident that the relevant momentum 
fraction is x and we are therefore lead to choose the physical variables to be given 

by 

^ ^(1 -Z)-K pI 

u = , X = — , K = — . (3.30) 

Indeed, the Lorentz transformation 

X k'jk^. , u{l-x) K^^K, , x k^^K, - K^^, 

(1 -U){X- U) Pa-Pj X-U Pa-Pj X-U Pa ■ Pj 

with K = Pa + Pj relates the above parametrization to one preserving the direction 
of the incoming parton, 

Rqa = -Pa (3.32) 

X 

1 — X 

Rq = upj + (1 — u) Pa — k± (3.33) 

X 

1 — X 

Rqj = (1 - u)pj + u Pa + k^ . (3.34) 

In order to derive the phase space convolution properties associated with the 
parametrization given above, we employ the formalism outlined in the appendix. 
Substituting 

u = r , X = ]- (3.35) 

w-\-y[l — w) w-\-y[l — w) 

the parametrization above is mapped to 



{-qa) ^ w{-pa) + {I - w) y Pj - qi_ (3.36) 
q^ {l-w){-pa)+w y pj + qx. (3.37) 

?i = (i-y)p.-, (3-38) 



with qj_ = —Sajyw{l — w) such that the generahzed phase space measure factors as 

dMlj, -QalQ) = d(/)2(pj, -Pa\Q) ■ (3.39) 

Having identified x to be the relevant momentum fraction we consider the 
hadronic colhsion in a frame where^ 

Pa-qa^-Pa-Pa , N ■ ^ -N ■ , (3.40) 

X X 

such that the phase space convolution properties of both parametrizations become 
equivalent at hadron level, 

d(l)2{(li:Q\Q;Pa:qa:Xa) = d^f (pj | Q; P„, p„, XX„) — ^^^^ _ - ^ ^^^ ' ^^'^^^ 

We stress that the crucial difference is related to the fact that, considering the 
physical forward evolution, our parametrization does generate a finite transverse 
momentum for the parton entering the hard process after additional parton emission. 
Averaging over azimuth, the initial-final splitting kernels take the form 

such that the splitting probability is 

dP - ^^{V{p^,z)} ^^^^^^ 6[x xa) _ _ ^ ■ (3.43j 

Note that in the coUinear limit, k, — > we have x — > z such that the coUinear 
behaviour is properly reproduced."^ 

The phase space boundaries are given by the requirement that Xa < x, 



1^ <P±< -, Z±^-\1 + Xa±{l-Xa)\ l 



21 ^ 'V (i-^a)s. 



(3.44) 

3.2.2 Initial State Spectator 

Initial state radiation with initial state spectator, Pa,Pb Qa,(l,Qb is described by 
the parametrization 

1 i^l — V — X 1 
Qa = — ■ — Pa H : Pb H : — ki_ (3.45) 

V + X X V + X V + X 

1 — V — X V 1 1, 

q = ■ Pa H ■ — Pb H ■ — k± (3.46) 

V + X X V + X V + X 

q,^ {'i- + l)pb, (3.47) 



^Note that there is no a priori relation between incoming hadron and parton momenta in our 
formulation. 

^For readability we have suppressed indexing a possible flavour change of the incoming parton. 
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Figure 2: Available phase space for a final-initial dipole with invariant momentum trans- 
fer ^Saj = ^/—i = 100 GeV and an infrared cutoff of 5 GeV. The shaded region is 
accessible starting at the hard scale, the region enclosed by the solid line is an example of 
the phase space excluded when starting at a lower scale. The phase space regions for an 
initial-final dipole are identical. For a final-initial dipole, the variables 



for the initial-final 
correspond to a collinear limit. 



1 — u. Note that in the latter case u ^ 1 and u ^ 



preserving q — Qa — qb = —Pa —Pb- The transverse momentum is defined to be purely 
spacelike in the dipole's rest frame and obeys 



kj_ = -(1 - V - X)-Sab 
X 



Sab = '^Pa ■ Pb 



The variables x and v are those introduced in [59], 

Pa-Pb Qa-q 



X 

qa- qb ' qa- qb 

and we define the physical variables to be given by 

z{\ — Z) — K K 

X = 



1-Z 

Note that the Lorentz transformation 



1-z 



Sab 



(3.48) 



(3.49) 



(3.50) 



Pb -Pa 



Pb -qa qa-Pa 



Pb ■ qa 



Pb -Pa qa- Pa 



-PaPa,u - 



qa-Pa 



{q^Pa. + P^alau) (3.51) 



does transform this parametrization to a parametrization where 

Sqa = —^Pa , Sqb = ^-^Pb ■ (3.52) 

X + V X 

Following the arguments of the previous section we then find the phase space con- 
volution 

d0f Pa-, qa, Xa, Pb, Qb, Xb) = 

( \Q;Pa,Pa,{^ + 't^)^a;Pb,Pb, —^Xb] ^ ^ „\ dp\ . (3.53) 

Averaging over azimuth, the initial-final splitting kernels take the form 
such that the splitting probability is 

with 

the ratio of incoming parton flux. ^ Note that in the coUinear limit, v, k — > and 
X — > z such that wc find the correct coUincar behaviour. 

We remark that it would be possible to keep the spectator unchanged upon 
properly substituting the integrations over the incoming momentum fractions (the 
Jacobian being equal to one). This, however, would invalidate the fact that the above 
given parametrization of the splitting kinematics does preserve energy-momentum 
locally involving the emitter-emission-spectator system only (in fact, after applying 
the relevant Lorentz transformation S this would constitute the inversion of the 
kinematics as used in the dipole subtraction context). A further argument to not 
keeping the spectator unchanged is that, following the discussion on soft and coUinear 
factorization, we see no reason why the emission off the colour connected system pa, Pb 
should leave ph unchanged except for a strictly soft and/or coUinear emission. 

The phase space limits are now determined from x > XaXb = r to be given by 




.± = ^|1 + x±(1-x),/1-^jJh1_) , p^<(L4)!M. (3.57) 



^As for the final state spectator, we have suppressed indexing a possible flavour change of the 
incoming parton. 
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Figure 3: Available phase space for emissions off an initial-initial dipole of mass 100 GeV 
with r = 0.02 and infrared cutoff 5 GeV. The shaded region is the available phase space 
when starting from the hard scale, the region enclosed by the solid line is an example of 
the phase space excluded when starting at a lower scale. 

4. Conclusions 

We have argued that parton showers, based on Catani-Seymour dipoles have a 
number of properties that turn out to be useful when one considers the matching with 
NLO matrix elements. Several other groups [57, 58] have already written a parton 
shower program, following this motivation. In these approaches, however, the choices 
of phase space boundaries and evolution variable were made rather intuitively. The 
question for coherence properties, in particular, was left open. 

Taking this as a starting point we have investigated the soft coherence properties 
of a CS parton shower formalism. We are lead to the transverse momentum of each 
dipole as the evolution variable of our cascade. Furthermore, we chose a hard scale 
that allows us to access the whole kinematically allowed phase space. We then 
explicitly show that such a parton shower can reproduce the expected Sudakov 
anomalous dimensions and hence include soft coherence effects properly. 

We have specified all details of such a parton shower that will be important for 
its implementation. In particular, we addressed the issue of transverse momentum 
from initial state radiation that will build up from several emissions in our case. An 
implementation of such a parton shower is underway. 
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A. Treatment of CoUinear Factorization 

In view of our probabilistic treatment of parton showers [50], we choose to rephrase 
coUinear factorization in a way that fits to the picture of assigning the usual phase 
space measure 

to each parton in the evolution, and not only final state partons. In turn, this 
allows us to derive phase space convolution properties using the general phase space 
factorization inherent to the parametrization of final state splittings with final state 
spectator. 

We therefore rewrite 

j dXa = J dXa J ^^^^^^ 6FiPa,Pa,Xa) (A.2) 



where 



SF{Pa,Pa,Xa) = -^5 ( ^) 6 ( ^-x. ) 9{x^{l - x^)) . (A.3) 



2-Pa ■ Pb \Pa' PbJ \Pa ' Pi 

Here, Pa denotes the lightlike momentum of the incoming hadron and the lightlike 
vector Pb defines the coUinear direction, i. e. is taken to be the momentum of the sec- 
ond incoming particle. We note that this extends straightforwardly to dimensional 
regularization. 

The phase space measure with an incoming parton then takes the form 

#n(Pl' ■■■:Pn\Q; Pa:Pa: ^a) = d0„+i(pi, ...,p„, -Pa\Q)SF{Pa: Pa: Xa) (A.4) 

in terms of the general phase space measure 

d0„+i(gi, = d0i(gi) • • •d0i(g„)(27r)^5(gi H h g„ - Q) ■ (A.5) 

The phase space convolutions as given in [59] take the form 

/ dXaf{Xa) dz d(l)n{pi,...,Pn\Q + ZQa) qa-PiJ{z)S{z - X)dx , (A.6) 

Jo Jo 



where qa = {l/z)pa = XaPa is constrained by coUinear factorization. It is then easy 
to prove that the same convolution within our framework reads 

dXafiXa) J —Pa-PiJ{z)6{z-x)d(j)^{pi,...,Pn\Q;Pa,Pa,ZXa)dx (A. 7) 



and gives rise to the factorization 

dXafiXa 



'^^^''^^\a ■ PiJ{x)6{x - Xa)^ ] d(j)^{pi, ...,Pn\Q; Pa,Pa,Xa) ■ (A. 8) 



X^ 
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